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NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


TECHNICAL NOTE D-798 

SOME SIMPLE SOLUTIONS TO THE PROB LEM OF 
PREDICTING BOUNDARY-LAYER SELF-INDUCED PRESSURES 
By Mitchel H. Bertram and Thomas A. Blackstock 


SUMMARY 


Simplified theoretical approaches are shown, "based on hypersonic 
similarity boundary- layer theory, which allow reasonably accurate esti- 
mates to be made of the surface pressures on plates on which viscous 
effects are important. The consideration of viscous effects includes 
the cases where curved surfaces, stream pressure gradients, and leading- 
edge bluntness are important factors. 


INTRODUCTION 


As the Mach number increases into the hypersonic range large pres- 
sure and temperature gradients can be induced on a flat plate by 
leading-edge bluntness and boundary- layer displacement effects. Any 
curvature to the plate surface simply adds still another complicating 
factor to the flow picture, and there may be any combination of these 
factors. In addition to the effect of these induced pressures on 
loading, the increased pressures can have an important effect on the 
skin friction and heat transfer on the plate surface . 

The purpose of this paper is to show the simplified theoretical 
approaches which are presently available and which allow reasonably 
accurate estimates to be made of the surface pressures on plates on 
which viscous effects are important and where surface curvature, stream 
pressure gradients, and leading-edge bluntness can be contributory 
factors . 


SYMBOLS 


a constant 


C “ li W T oo/^oo T W 



diameter of cylinder 
f(T w ,T 

r>y > H») in equation for laminar boundary- layer growth 
(see eqs. ( 3 ) and (l4)) 

index number (0 for planar model, 1 for cylindrical model) 
hypersonic similarity deflection angle (M^ • deflection angle) 

hypersonic similarity deflection angle using slope at edge of 
boundary layer 

coefficient in hypersonic- similarity- theory- laminar-boundary- 
layer growth equation (see eq. (2) and fig. l) 

nose drag coefficient 

Mach number 

exponent in power equation for surface coordinates of curved 
plate ^y « x m ^) 

Knudsen number 
Prandtl number 

exponent in power equation for surface pressure variation 
with x 

ratio of local surface pressure to free- stream static pressure 

pressure 

Reynolds number 

Reynolds number based on undisturbed free- stream conditions 
and leading-edge thickness or rod diameter, respectively 

absolute temperature 

leading- edge thickness 

distance along axis of flat plate from leading edge 
model span 
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a 

7 

5 

e 

IJ- 

I 

X 


angle of inclination of flat plate relative to free- stream, 
flow direction 

ratio of specific heats 

"boundary- layer thickness (see eg.. (2)) 

correction factor for y in hluntness-induced- pressure 
equation (see eq. (^l)) 

dynamic viscosity 

coefficient in induced-pressure equation for rod (see eq. ( 50 )) 


viscous interaction parameter 



as local surface slope 

Subscripts: 

b nose bluntness induced 

1 inviscid 

r adiabatic wall 

v viscous induced 

w wall 

x surf ace di stance 

oo free stream 

0 stagnation 

2 local inviscid pressure 

A prime designates the first derivative with respect to the subscript. 

Superscripts in parentheses refer to the order of an equation devel- 
oped in series. 
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SOME SPECIAL SOLUTIONS FOR TEE SELF- INDUCED 
PRESSURE GRADIENT 


Useful results may be obtained if the similarity theory for a 
boundary layer in hypersonic flow (ref. l) is used as a basis for 
predicting self-induced pressures to a higher order than may be obtained 
from the so-called "weak- interaction" or "strong- interaction" solutions. 
(See ref. 2 for a discussion of these theories.) First, consider the 
boundary- layer thickness itself in hypersonic flow. To satisfy the 
requirements of this hypersonic similarity theory, a pressure variation 
at the edge of the boundary layer is postulated as follows : 


3?W 


« x 


n 


( 1 ) 


Then according to the hypersonic similarity theory (ref. l) the boundary- 
layer thickness is (assuming Prandtl number unity in the derivation) 


5 


GK 4 ( 7 ,n) 


Mro /C 


. ; P 



( 2 ) 


where 


G = 1.7208 y -- ~ 


tT °‘ 5859 


(3) 


and K 4 (7,n) is the coefficient determined from hypersonic similarity 
theory. The variation of K4 with negative exponents in the pressure- 
variation law (favorable pressure gradients) is shown in figure 1. The 
curves shown are different from those in reference 1 because additional 
values from reference 3 of the integral function used to obtain K 4 have 
been utilized to obtain more accurate fairings of the integral function. 

If hypersonic tangent-wedge theory (ref. 4) with the slope of the 
outer edge of the boundary layer as the effective wedge angle is assumed 
to be sufficient to determine local surface pressures then 


'00 


p = fK 8 i) - f (%) 
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The boundary- layer slope is upon differentiation of equation (2) (with 
Ki|. assumed to be constant) 



Jc 

2 J^co } X 


x dP 
P dx 


(5) 


or in terms of the Lees-Probstein interaction parameter X (ref . 5) 
where 



(6) 


the hypersonic similarity deflection angle 


may be written 


K 

2 




i + 


2P 



(7) 


Wow consider which values of the pressure-law exponent n are of 
interest if self- induced pressures only are to be considered. The 
variation of the local value of the exponent n with self- induced pres- 
sure ratio P is shown in figure 2. The value of n is zero at P = 1 
and at increasingly large pressures the curve of n is asymptotic to 
-0.5, the value obtained from strong- interaction theory. The values of 
n shown in_figure 2 were obtained by assuming that the pressure is 
linear in X with the following form (developed in more detail in a 
later section) 


P = 1 + aX 


(8) 


since for any function (of x, say) the local value of 
power- law fit to the function may be expressed as n = 

equation (8) f(x) is (l + aX) and thus the exponent 
of pressure ratio is given by 


the exponent in a 
than for 

f(x) 

n as a function 


1 P - 1 

2 P 


n 


(9) 
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in helium (7 = 5 / 5 ) the coefficient K 4 is 

generally less than 2 percent) in the range 
(7 = 7 / 5 ) the maximum deviation from unity is 

Weak Interactions 

First consider equation ( 7 ) in the light of weak- interaction 
theory. From hypersonic small-perturbation theory (see ref. 4) 


Figure 1 shows that 

essentially unity (error 
0 > n > - 0.5 and in air 
about 10 percent . 


p( 2 ) = 1 + 7 K & + 2lZ_t_^K§ (10) 

Assume that a first-order correction for the local pressure is sufficient 
in equation ( 7 ) • That is, the value of P to be inserted in equation (10) 
will be obtained from the first-order solution to the weak- interaction 
boundary- layer problem (ref. 5 ) where 

pt 1 ) = 1 + Z®^2. (11) 

Substituting equation (ll) into equation ( 7 ) with = lj and then this 
combination into equation (10) results in the following relation 



( 12 ) 


Equation (12) may be compared with the original Lee s-Prob stein equa- 
tion (ref. 5) for second-order weak interaction which in the present 
nomenclature may be written as follows: 
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P 


( 2 ) = 1 + & 



(15) 


Equation (13) may tie derived readily from equation (lO) "by utilizing 
equation ( 7 ) and n = 0. 


Figure 3 has teen prepared to show the difference between the use 
of equations (12) and ( 13 ) in the case mentioned previously where the 
boundary- layer thickness coefficient is essentially unity, that is, 

the flat plate in helium. In addition to the Prandtl number one solution, 
curves for Prandtl number O .725 are also shown for which the preceding 
development is taken to hold except that G is given by the following 
equation adapted from reference 6: 


G = 1.648 ^- 1 |^ + 0.352 


\T r 


(14) 


The reduction in induced effects due to the first-order correction 
for local pressure is certainly significant and appears to be corrobo- 
rated by experimental results obtained by Henderson and_Johnston on a 
5° wedge tested in helium (ref. j ) . In figure 3 both X and the 
pres sure- difference ratio are based on the inviscid conditions that 
would occur on a wedge having a sharp leading edge were there no boundary 
layer present. For the two stations nearest the leading edge the experi- 
mental data are shown both as measured and with an increment in pressure 
subtracted to account for leading- edge bluntness . This increment 
ascribed to leading- edge-bluntness effects is estimated from a correla- 
tion of characteristics solutions for the flat plate with a sonic-wedge 
leading edge at angle of attack (ref. 8) by parameters suggested by 
Chernyi (ref. 9) and based on the assumption that the boundary- layer- 
induced and bluntne s s- induced effects may be treated independently. The 
effect of leading-edge bluntness is indicated to be small. The 10° wedge 
data also shown on this plot were obtained at too low a value of X to 
aid in distinguishing between theories considering data uncertainties. 


Strong Interactions 

Thus far the discussion has been restricted to considering rela- 
tively small induced pressures, that is, restricted more or less to the 
weak- interaction regime of self-induced pressures. Assuming that the 
tangent-wedge approach is still applicable, a more general approach, 
insofar as the basic pressure equation used is concerned, will now" be 
attempted. In this case the hypersonic equation given by Lees (ref. 7) 
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will be used. This equation is in terms of boundary- layer slope 


( 15 ) 


Based on equation (7), with Kij. = 1 and with equation (15) , calcula- 
tions have been carried out for helium. These calculations were made by 
iteration for a given value of X by assuming for each iteration that 
the pressure ratio is a linear function of X starting with the result 
from strong- interaction theory. (See refs, 4 and 1.) In this way the 
results are easily calculated. Iterations are obtained by averaging the 
assumed pressure ratio and the pressure ratio resulting from this assump- 
tion. This average is used to start the next iteration. The results are 
shown in figure 4. The convergence is rapid. The use of the more accu- 
rate pressure law (eq. (15)) makes only a small difference from the 
results obtained by using the expansion of equation (15) for large Kg or 



The differential equation which results when equation (l6) is substituted 
into equation (7) may be solved to give the asymptotic solution known as 
strong- interaction theory 



where G is given by equation (3) and the required values of K4. are 
given in figure 1. A comparison of the complete hypersonic equation 
(eq. (15)) with its "strong" shock approximation (eq. (l6)) and its 
"weak" shock approximation (eq. (10)), is shown in figure 5. 

The result of the iteration procedure utilizing equation (15) and 
equation (7) will hereafter be referred to as the "complete theory" (for 
the sake of briefness even though fundamental completeness is not 
implied). The results from the complete theory are the same as those 
from equation (12) for values of GX less than approximately 3. For 
GX greater than approximately 1 the curve for the variation of P with 
GX is essentially linear and may be accurately represented by the 
equation 


P = 1 + 


7(7 + b) v 2 


•5 + ? K 6 


7 + 1 

T 4 
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P 


= O.85 + 


3 fim) G x 
k v 2 


(18) 


Interestingly enough, for practical purposes this result could have been 
obtained from equation (17) if K^ were assumed to be unity and 

p(°) = P - 1, even though equation (l6) deviates markedly from equa- 
tion (15) at the lower pressure ratios. For P > 4 the induced-presssure 
ratio can in general be taken as 


P « 1 + 


3 h(v + 1) 
k\J 2 



n=-l/ 2 


(19) 


while for P < 4 equation (12) should be adequate and, in general, data 
should lie in the region defined by equations (12) and (19) • A compari- 
son between equations (12), (13) , (19) , and the complete theory (essen- 
tially eq. (l8) is shown in figure 4. 

Consider now the insulated plate in air where the value of the 
coefficient deviates by an appreciable amount from the value unity. 

(See fig. 1.) Here the fact that is greater than 1 offsets, to 

some extent, the result of P > 1 on the boundary- layer thickness. (See 
eq. (2).) Thus in the low X range the actual induced pressures will be 
expected to be between the values given by equations (12) and (13), tat 
closer to the values from equation (12) . As X increases the actual 
pressure would be expected to approach the complete theory with a con- 
stant value of Klj. = 1.1 which is essentially equation (19) • 

These various theories are shown in figure 6 with applicable data. 
Shown are Kendall's data from an insulated flat plate at M ra = 5*8 
(ref. 10) and data obtained in the Langley 11- inch hypersonic tunnel at 
= 9.6 on a noninsulated plate with a temperature gradient for which 
the data have been corrected to the insulated plate case as described in 
reference 2. At ^ = 9.6 data both with the plate surface alined with 
the flow and at an angle of attack have been obtained. (A portion of 
this data has previously been reported in refs. 11, 12, and 13 •) The 
data for zero angle of attack extend the data given in reference 2 and 
are believed to be more accurate while angle of attack has been used to 
extend the range of X and to give an additional check on the validity 
of the correction to insulated plate conditions. At angle of attack the 
local inviscid Mach number departs appreciably from 9*6 as shown i-n 
figure 6 (b) . The original data were picked at a time when the temperature 
distribution on the plate was the same for all Reynolds numbers and 
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angles of attack. This temperature distribution is shown in references 11 
and 12 . 

At zero angle of attack (fig. 6 (a)) there is good agreement between 
Kendall's data and present data in the range where they overlap. The data 
points in general are slightly above the curve for weak- interaction theory 
with the local pressure correction (eq. ( 12 )) but approach the curve for 
"complete theory" with = 1.1 at the higher values of X. The pressure 
data at angles of attack (fig. 6 (b)) show good_agreement with the data at 
zero angle of attack, at least for values of X up to about 10 , and indi- 
cate that the complete theory with = 1.1 (dr eq. (19)) is adequate at 

least in the range of X shown. 

BOUNDARY- LAYER- INDUCED EFFECTS IN A NON- SELF- INDUCED 

PRESSURE GRADIENT 


Strong Non- Self- Induced Pressure Gradient 

Effect of wall curvature .- First, with equation (l) as a requirement, 
the slope of the edge of the boundary layer (eq. (k ) ) in hypersonic simi- 
larity form may be written as 


K& = = § T~~ — ( 1 - n)K k 

/P2/P00 


( 20 ) 


Then assume that the boundary layer is growing on a surface of sufficient 
curvature so that equation (l6) well represents the pressures induced by 
the surface slope alone. In general, a change in the slope of the sur- 
face induces a pressure rise above the pressure that would have existed 
as follows (from a finite perturbation on eq. (16)): 


P w- P 2 


= 2 (ak)k + (ak) : 


(21) 


where K = if £K = K§ and K5 is assumed to be small then equa- 

tion (21) with equations (20) and (l6) gives 


ro voi vjj H f 
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P w~ P 2 


2(7 + 1) 


GXjl - n)K^ 


( 22 ) 


which, since 2(AK)K » (£K) P , holds for 


£ » JHr - n)K ^ 


where n is obtained from 



in the form of equation (l6) . 


Equation (22) does not reduce to equation (17) when n = -l/2, since 
the self- induced effect is derived as a small perturbation on a pre- 
existing pressure gradient. 


If the ratio of the local wall pressure, with the self-induced effect 
included, to the inviscid local pressure is desired equation (22) may be 
written as 


— = 1 + f ~ — fM 1 - n)K k (25) 

p 2 V 7 2 p 2 / p oo 


which has the same limitation as equation (22) . 


Now for hypersonic flow 


£2 = EX 

V T 2 P 2 


(24) 


For an oblique shock compression with pressure ratios not too close to 
unity 


£2 = 7-1 £2 + 
7 + 1 Poo 


1 


(25) 


If, further, the unity term in equation (25) can be neglected then 


£2 _ ITT T 

\o " p 2V7 - 1 


(26) 



12 


Then in terms of inviscid local conditions equation (23) with equa- 
tion (26) becomes 



( 27 ) 


Equation (23) may also be presented in terms of surface slope by using 
equation (l 6 ) or 


-w 


= 1 + 


/ 2 

\jv(7 + 1) K2 


•(1 - n)K 4 


K 2 » ft. 


7(7 + 1 ) 


■(1 - n)Ku 


(28) 


The externally imposed pressure gradient .- The case where the pres- 
sure gradient is externally imposed may also be simply treated. Again it 
is assumed that the increase in surface pressure is small compared to the 
inviscid local pressure. The local pressure ratio is assumed to be given 
by the first-order form of equation ( 10 ), which is 


r = 1 + 7K& > 2 

where the required form of equation (7) is 


(29) 


%,2 


M 2 6 i - f 


M, 


/V 


i? 1 ' 


(30) 


and with isentropic external flow or isentropic compression on a curved 
plate 



X 

00 

27-I 


(1 



n)K 4 


(31) 


rov>j v* (-> c 
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in which the variation of the local inviscid pressure ratio pp/p^* with 
surface distance has the form of equation (l) and the value of (l - n)K k 

is obtained from this pressure law. In this case a> refers to some 
convenient reference condition such as the local free- stream Mach number 
at the leading edge. Equation (29) with equation (31) gives the result 



1 + f Vi (1 ' n)Kl > 

Pg\2? 


(52) 


With reference to equation (32) the following point is made; that is, 
by hypersonic approximation for an isentropic compression, equation (2^) 
is 


X„ 


27-I 




(33) 


so that equation (28) may be written as 


( 1 ) 


r w 


P 2 


= 1 + 


2|x 2 (i - n ) K 4 


( 3*0 


Equation ($k) has the same form as equation (27) and gives not too dif- 
ferent a result even though the development of the two equations was 
somewhat different. Note that, with n = 0 and p 2 = equation (3^) 

reduces to equation (13) (even to second order, it can be shown) . Thus 
equation (3^) is essentially the form of weak- interaction theory for 
self-induced effects of a boundary layer growing in a non- self-induced 
pressure gradient and the term (l - n)Kij. is the correction term due to 
the external gradient . 


1 
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n 

(l - 

n)K^ 

7 = 7/5 

7 = 5/3 

0 0 

1.000 

1.000 

-.1 

1.119 

1.104 

- .2 

1.237 

1.208 

-•3 

1.364 

1.313 

-A 

1.502 

1.421 

-.5 

1.649 

1.533 

- .6 

1.801 

1.648 

-.7 

1.965 

1.769 

-.8 

2.145 

1.895 


Shown here are values of 
(l - n)K^ for some arbitrary 

values of n for favorable 
pressure gradients on an insu- 
lated plate in both air and 
helium. Other temperature 
ratios and values of n for 
air and helium or other values 
of 7 may easily be calcu- 
lated by using figure 1 and 
reference 1. 




Weak Non- Self- Induced Pressure Gradient 

It is clear that with a power-law variation of pressure with sur- 
face distance the boundary- layer growth and hence the self-induced 
effects depend upon the local conditions with a generally small, but 
not necessarily negligible, correction due to the exponent in this 
power- law variation of pressure. This dependence on local conditions 
will now be used to investigate the case where the change in local 
pressure, whether by external causes or changes in surface slope, is 
small compared to the local self-induced pressure . Based on inviscid 
local conditions, with the exponent in the pressure law having negli- 
gible deviation from n = -l/2 (the strong- interaction case), equa- 
tion (17) may be written 


fw 

P 2 


3 

4 


7 ^ 7 + ^ K4GX2 


(35) 


where the subscript 2 indicates local inviscid conditions and the sub- 
script w indicates surface conditions. Equation (35) written in terms 
of reference conditions designated by the subscript <» is 



7(7 + 1 ) 
2 



( 36 ) 


With hypersonic flow and isentropic pressure changes, equation (33) may 
be applied from which 


H KMCW 




2E 


1 


^2 £2 = fei \ 27 

*oo P°° \ °°/ 


15 


( 57 ) 


and for small pressure deviations from the reference conditions 


^2 ^2 ! + 1_ Al? 

V P 27 P 

Aqq 00 ' CO 


(38) 


where Ap = - p and, to the present approximation, equation ( 38 ) 

applies equally well for pressure changes caused by weak shocks. Equa- 
tion ( 38 ) with equation ( 38 ) becomes 


1 + 1 . 

2 ? p <* 7 - 1 K, 

in which p^^ designates the strong- interact! on self- induced pressure 

calculated for free- stream conditions which are the same as the refer- 
ence conditions 00. if the inviscid pressure change Ap^p is due to 

a small surface inclination, then in terms of the hypersonic similarity 
parameter K = Mjn 



Pw 

P2 



(HO) 


Equations (39) and (40) allow an estimate to be made of such factors as 
(l) the effects of small errors in evaluating free- stream Mach number 
or in setting the angle of attack of a flat plate, ( 2 ) the effects of 
small body curvatures, and ( 3 ) the possible effects of external pressure 
gradients. Factors (l) and ( 2 ) are straightforwardly obtained from 
equations (39) and (4-0) . Factor ( 3 ) is interesting because of the pos- 
sibility of using equation (39) to evaluate the effects of longitudinal 
pressure gradients in conical nozzles such as those which have been used 
to provide much of the information on self-induced effects in helium. 

Data have been obtained in 1-inch- and 2- inch- diameter conical helium 
nozzles at Langley Research Center (refs. 7 and l4) for which the Mach 
number increases about 0.5 in the length of the measuring stations on a 
flat plate. These tests were performed at Mach numbers between 16 and 24. 
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In the case of Erickson's tests (ref. 14) the correlating parameters were 
evaluated at tunnel-empty local conditions, corresponding to the sub- 
script 2 in the preceding development. In the investigation by Henderson 
and Johnston (ref. 7) "the correlating parameters were evaluated at the 
tunnel-empty conditions at the leading edge of the model. Equation (39) 
indicates a difference between the two methods of evaluation of the slope 
of P = P(X) of at most about 4 percent. Also, the P(x) line from a 
leading-edge-condition evaluation will not pass through the P = 1, 

X = 0 point. 


THE BUJETNESS- INDUCED— VISCOUS- INDUCED INTERACTION PROBLEM 


Theoretical Approach 


The problem of the interaction between bluntness- induced and 
boundary- layer- induced flow fields on a flat plate is of considerable 
importance. One approach to this problem is suggested by the work of 
Baradell and Bertram (ref. 8). In this paper, by means of the charac- 
teristics theory, surface pressure distributions on curved plates with 
surface coordinates given by the equation y « x m were investigated. 

It was found that addition of the pressure increment due to local plate 
angle (tangent- wedge theory) and the pressure increment due to the 
blunt leading edge (as on a flat plate) gave a result which agreed well 
with the characteristics theory. This blunt-leading-edge pressure-ratio 
contribution for p^/p^ = 1.2 is given accurately by 


*b 


= 0.l87e 


y/7(r - i) 


M^k 


1 2 /3 


x/t 


+ 0.74 


m 


where e is a small correction term to the function of 7 shown above 
which can be represented by e « 1 - [b. 0048/(7 - l)^] • Equation (^l) is 
found to give good re stilts from constant 7 of 6/5 to 9/5 based on the 
results in reference 8. 

The curved plate in the problem of reference 8 is replaced by the 
displacement thickness of the boundary layer. The growth of this 
boundary layer is assumed to be controlled by the pressures on the plate 
and by the undisturbed free-stream Mach number; that is, no account is 
taken of the entropy jump at the leading edge . The decrease in Mach num- 
ber associated with the leading-edge entropy jump is assumed to be com- 
pensated for by the decrease in Reynolds number. 

Two different cases for the problem of surface pressures induced by 
blunt- leading- edge viscous interaction will be shown. In the first case 
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the blunt-leading-edge-induced pressures are assumed to be large com- 
pared -with the viscous- induced increment in pressure. In the second 
case the viscous-induced increment in pressure is too large to be ignored 
and contributes significantly to the induced effect on boundary- layer 
growth. If the blunt- leading-edge- induced pressures Eire assumed to be 
large compared with free- stream pressure and much larger than the viscous- 
induced increment, then the slope of the edge of the boundary layer in 
hypersonic similarity form may be given by equation (20) in the following 
form (n = -2/3 from eq. (4l)): 


k 5 = KP* = 


6 /Tip 


(42) 


Taking the coefficient as unity (see fig. l) and substituting equa- 

tion (42) into the first-order form of equation (10) to represent the 
tangent-wedge portion of the solution yields 


Ap^ 1 ^ ^ 5 7®^ 

p “ 


( 43 ) 


or to the second order 


Ap 


( 2 ) 


v 


Ap, 


( 1 ) 


P„ 


7 + 1 /Ap, 

4y 


( 1 )' 


(44) 


The parameter that controls the blunt-leading- edge- induced pressure 
in equation (hi) may be written as follows: 


I^k = kR t -2 


( 45 ) 


Thus for a given value kR^I^C the wall pressure (p w = p^ + Ap v ) is 

a function of only. This result can be obtained from reference 15 

which uses a considerably different development. 
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In the second case, the effect of the viscous- induced pressure 
increment on the boundary- layer growth must be accounted for. In this 
case, equation (20) is written as 


k 5 


(1 - n)K 4 
2 


GX„ 




(46) 


The form of equation 



(15) for large values of K5 is taken as 
and when this equation is substituted in equa- 


tion (46) the solution is 


Ap 

p. 


z = 1 < 


K\ + ZhjLii[(l - n)K^] 2 (a>cJ 2 - 
V V 00 / 


Pb 


(47) 


Comparison With Experiment 

In order to check these predictions of viscid- inviscid interaction 
effects figures 7 and 8 have been prepared which show much of the data 
available on flat plates in air and helium. In figure 7 data obtained on 
essentially insulated flat- leading- edge plates in helium from references 7, 
14, 16, and 17 plus some new data obtained in the Langley 11- inch hyper- 
sonic tunnel (see appendix) have been put into groups with approximately 
the same value of kR^M^C. First it may be noted that there appears, in 
general, to be good correlation of the data for given values of the param- 
eter kRt/M^C . Some random deviations from correlation may be noted} how- 
ever, these deviations do not appear to be associated with any particular 
Mach number or Reynolds number trend. The sort of deviation or scatter 
referred to is illustrated in figures 7(^) and 7( c ) "by the data obtained by 
Henderson and Johnston (ref. 7) for values of kR-t/^C of about 5 (square 
symbols) and 6.6 (diamond symbols). These sets of data (compare flagged 
and unflagged symbols) were obtained under almost the same conditions. 

The strong viscid- inviscid interaction theory (eq. (47)) does a good 
job of predicting the surface pressures except perhaps at values of 
kR t h&p below about 1.5 where the pure viscous theory appears to be 
preferable. At values of kRtM»C where inviscid effects are dominant 
(above about l4) the first-order weak viscid- inviscid interaction theory 
(eq. (43)) gives much the same answer as equation (4l) and both agree 
well with the data. 



19 


Much the same result is obtained for the data available in air 
flows shown in figure 8 . The data (refs. 17 to 22 plus new data 
obtained in the Langley 11- inch hypersonic tunnel), in general, show 
good agreement with the theories. It should be noted that, although for 

the data presented in figure 7 the lowest value of kR^/M^C is about L, 
the data presented in figure 6 for comparison with pure viscous theory 
had very much lower values of kR^M^C, values which varied from 0.15 

to 0.5 for the M = 5.8 data (ref. 10) and from 0.05 to 0.1 for the 
M = 9-6 data of the present tests. 

There are, however, some cases where there is a major disagreement 
between experiment and theory. This disagreement generally occurs for 
data obtained in the region close to the leading edge. Data obtained in 
approx ima tely the first two diameters from the leading edge in figure 8 
are indicated by the filled symbols which generally correspond to the 
data showing the large disagreement. It is in this region close to the 
leading edge that the effect of nose shape is important in determining 
the pressure distribution. 

Analysis of data in the near- leading- edge region, however, can 
also be complicated by low-density effects. This is illustrated in 
figure 9 by data obtained in the University of California low density 
tunnel by Aroesty (ref. 25) at M ~ 3*9 and by Schaaf, Hurlbut, and 
Talbot (ref. 18) at M « 5.7. Aroesty' s data (fig. 9(a)) were obtained 
at very low Knudsen numbers for the leading- edge thickness and_shows 
qualitatively the same slope pattern at the larger values of X as is 
shown in the pressures measured by Schaaf, Hurlbut, and Talbot (fig. 9("b) ) 
at Knudsen numbers which ranged from relatively high (2.1) to quite low 
(O.Ol) . Figure 10 can be used to illustrate that this rather similar 
behavior of the data from reference 18 for the widely different Knudsen 
numbers is fortuitous. In figure 10 the pressure data of reference 18 
are plotted against x/t rather than against the parameter X for the 
first 10 diameters downstream of the leading edge of the plate. Shown 
also for comparison purposes are Creager's essentially inviscid data 
(ref. 20) obtained at a relatively high density on a flat- leading- edge 
plate. Note that except for the magnitude the shape formed by the 
University of California data for the lower Knudsen numbers (N^ ^ O.OL) 
is much the same as that of Creager's. The University of California data 
at the higher Knudsen numbers (Nj^ of 0.1 and 0 . 2 ) is not data from the 
near- leading-edge region and only coincidentally has roughly the same 
level of pressure as the data from the first two or three diameters at 
the lower Knudsen numbers. This fortuity is further corroborated by the 
fact that even though Aroesty and Creager obtained their data at about 
the same Mach numbers the maximum pressure ratio obtained by Creager is 
about 2 (fig. 10) but by Aroesty is about 3*5 (fig* 9(a)). 
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Oguchi has considered the problem of the prediction of the viscous- 
induced pressures in this near- leading- edge region (ref. 24). By uti- 
lizing a wedge- like boundary- layer growth in this region, apparently he 
has obtained a solution which successfully predicts the pressure levels 
obtained at low densities in references 23 and 25* 

Extension to the Case of Blunt-Nosed Rods 

Equations (43 ) and (47) may be written in terms of the inviscid param- 
eters for both flat plates and rods if for the rod the boundary layer is 
not too thick compared with the radius of the rod and results are obtained 
far enough back of the nose so that the three dimensionality of the nose is 
not important. For example, if the nose of the rod has a hemispherical 
shape then the number of diameters downstream from the nose must be large 
compared with jt/8 (ref. l) . 

This more general form of equation (43) is 


and of equation (47) is 



in which j = 0 for the planar case and j = 1 for the rod alined with 
flow (see eq. (5) of ref, 8). The symbol t represents the leading-edge 
thickness in the two-dimensional case and the diameter of the rod in the 
cylindrical case. 

Calculations based on equation (49) are presented in figures 11 and 
12 for air and helium. Figure 11 shows results for the two-dimensional 
case utilizing equation (49) with equation (4l) . The solid line is the 
inviscid case given directly by equation (4l) . The increments between 
the solid curve and the dashed curves represent the viscous- induced" pres- 
sure increments predicted by equation (49) „ The parameter G (eq. (3) 
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or eq. (l4) ) takes into account the wall temperature level. In computa- 
tions for figure 11, the value of n in equation (49) was taken as con- 
stant at a value of 0.6. The error that this introduces can be assessed 
by noting that the local value of n is the local slope of the curves 
in figure 11. 


To generalize from the experimental insulated- plate results pre- 
sented previously and compared to equations (48) and (49) (or to their 
less general forms - equations (43) and (47)), consider values of the 

parameter G When the value of this parameter is greater than 

about 2/3, the effects of leading-edge bluntness are virtually negligible. 
When the value of this parameter is less than about 1/2, blunt- leading- 
edge- induced pressure effects are important and the theory developed in 
this paper gives a good prediction of the surface pressures. For values 
of the correlation parameter less than about l/lO the viscous-induced 
contribution is small to negligible compared with the contribution of 
the blunt leading edge to the surface pressures. 


An analogous presentation is made for the blunt rod in figure 12. 
In this case equation (49) was used with an empirical equation obtained 
by Van Hise (ref. 26) based on characteristics solutions for rods with 
various bluntnesses. This equation is 




0.55 


(50) 


in which | = 0.060 for 7=7/5 and £ = 0.075 for 7 = 5/5- The 
solid line in figure 12 is the inviscid case given by equation ( 50 ) • 
Again the increments between the solid curve and the dashed curves 
represent the viscous-induced pressure increments predicted by equa- 
tion (49) . In computations for figure 12, the value of n in equa- 
tion (49) was taken as that given by equation ( 50 ) . (See discussion 
before eq. (9)-) Th is introduces some error since, judging from the 
computed curves, the value of n for the case with viscous interaction 
will always be less than the value of n given by equation ( 50 ) . In 
figure 12, d has been used to represent the t in equation (49) as in 
equation ( 50 ) . 


CONCLUDING REMARKS 


This paper has shown the simplified theoretical approaches, based 
on hypersonic similarity boundary- layer theory, presently available 
which allow reasonably accurate estimates to be made of the surface 
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pressures on plates on which viscous effects are important. The con- 
sideration of viscous effects includes the cases where curved surfaces, 
stream pressure gradients, and leading-edge hluntness are important 
factors . 

In a simple attack on the problem of the effect on the surface' 
pressure of the interaction between flow fields induced by leading-edge 
blunting and boundary- layer displacement effects, equations have been 
derived which in general give a good prediction of the surface pressures. 
The parameter determining the level of the surface pressure induced on a 
two-dimensional flat plate by the interaction of blunt-leading-edge- 
ind uced and boundary- layer- induced flow fields was found to be 

where G is a wall temperature and Mach number function in 

the laminar-boundary- layer growth equation, k is the nose drag coeffi- 
cient, is the Reynolds number based on leading-edge thickness, 

is the free- stream Mach number, and C is the coefficient in the linear 
viscosity law. Correlations of the available data according to this 
parameter indicate the effects of leading-edge bluntness are virtually 
negligible when the value of this parameter is greater than about 2/3. 
When the value of this parameter is less than about l/2, blunt-leading- 
edge- induced pressure effects are important and the theory developed in 
this paper gives a good prediction of the surface pressures. For values 
of the correlation parameter less than about 1/10 the viscous- induced 
contribution is small to negligible compared with the contribution of 
the blunt leading edge to the surface pressure. 

An extension of the theory to the case of the blunt-nosed rod 
alined with the free- stream flow also allows a prediction of the surface 
pressure induced by viscid- inviscid interaction for this special three- 
dimensional case . 


Langley Research Center, 

Rational Aeronautics and Space Administration, 
Langley Field, Va,, February 9> 19^1 . 
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APPENDIX 

DESCRIPTION OF NOZZLES AND MODEL USED TO OBTAIN DATA IN THE LANGLEY 

11-INCH HYPERSONIC TUNNEL 


Tunnel and Nozzles 

The unpublished data shown in this report were obtained from tests 
conducted in the Langley 11-inch hypersonic tunnel. The Mach number of 
this blowdown facility can be varied by changing nozzles, in the present 
tests, nozzles giving nominal Mach numbers of 6.8 and 9*6 with air and 
18 with helium were used. The Mach number 6.8 air nozzle is a contoured 
two-dimensional nozzle machined from Invar to minimize deflection of the 
nozzle throat due to thermal gradients. A calibration of this nozzle 
may be found in reference 27. The Mach number 9*6 nozzle is a contoured 
three-dimensional nozzle with square throat and test section. The cali- 
bration for this nozzle is given in reference 2. The Mach number 18 
helium nozzle is contoured to obtain uniform flow and has a circular 
cross section. A complete description and calibration of this nozzle 
and other helium nozzles used at the Langley 11-inch hypersonic tunnel 
may be found in references 28 and 29. 


Models and Instrumentation 

The model used in this investigation is the slab wing with the 
various leading edges shown in figure 13. The 36° wedge leading edge 
was designed to give sonic velocity behind the shock at Mach number 18 
in helium and the ^3° wedge leading edge was designed to give approx- 
imately sonic velocity behind the shock at Mach numbers of 6.8 and 9*6 
in air. Surface pressures were measured by means of the aneroid-type 
six-cell recording units described in reference 30. The skin tempera- 
tures were obtained from ehromel-alumel thermocouples formed from No. 30 
wire. The thermocouple wires were spotwelded to the underside of the 
surface and the skin temperatures were recorded on self-balancing 
potentiometers. Some typical wall -temperature distributions are shown 
in figure l4. These temperature distributions are ones obtained with 
the 0.060 -inch -thick flat -face leading edge; however, they are represen- 
tative of the distributions obtained for all of the various leading 
edges. For the nominal Mach number 6.8 tests, temperature distributions 
are shown for two different Reynolds numbers. For the tests at Mach 
numbers of 9*6 and 18 the effect of Reynolds number on the temperature 
distribution was small and only one distribution for these Mach numbers 
is shown. 
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Figure 1.- Variation of coefficient in hypersonic similarity boundary- layer- thickness law 

with pressure gradient. 



Figure 2.- Variation of exponent n in power law 
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Figure 3.- Comparison of theory with experiment on an insulated flat plate in helium 
for weak-interaction theory and a higher order approximation. 





Figure 5.- Static-pressure ratio across oblique shock as a function of the 
similarity deflection angle according to hypersonic similarity theory. 
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Figure 6.- Concluded 
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(a) 


kR t 

— in range from about 4 to about 9. 


Figure 8.- Correlation of induced pressures on flat plate in air and 
comparison of theory with experiment. Solid symbols indicate data 
obtained two diameters or less from leading edge. 
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figure 9.- Effect of low-density flow on induced pressures in near-leading- edge region 

of flat plate . 
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Fig-ure 10.- Ratio of measured wall pressure to static pressure as a function of 
number of diameters downstream of flat leading edge of plate . 




figure 11.- Generalized chart shoving effect on ratio of vail pressure to static pressure 
of interaction of boundary- layer displacement and blunt- leading- edge- induced effects 
on flat plate . 
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Figure 12.- Concluded. 









